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Abstract. We have performed p-process simulations with the most recent stellar 
(n, 7) cross sections from the "Karlsruhe Astrophysical Database of Nucleosynthesis 
in Stars" project (version v0.2, http://nuclear-astrophysics.fzk.de/kadonis). The 
simulations were carried out with a parametrized supernova type II shock front model 
("7 process") of a 25 solar mass star and compared to recently published results. A 
decrease in the normalized overproduction factor could be attributed to lower cross 
sections of a significant fraction of seed nuclei located in the Bi and Pb region around 
the 7V=126 shell closure. 
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1. The u p processes" 

A a p process" was postulated to produce 35 stable but rare isotopes between 74 Se 
and 196 Hg on the proton-rich side of the valley of stability. Unlike the remaining 99% 
of the heavy elements beyond iron these isotopes cannot be created by (slow or rapid) 
neutron captures [1] , and their solar and isotopic abundances are 1-2 orders of magnitude 
lower than the respective s- and r-process nuclei [2j [3]. However, so far it seems to be 
impossible to reproduce the solar abundances of all p isotopes by one single process. In 
current understanding several (independently operating) processes seem to contribute. 

The largest fraction of p isotopes is created in the "7 process" by sequences of 
photodissociations and (3 + decays [U [5j [6]. This occurs in explosive O/Ne burning 
during SNII explosions and reproduces the solar abundances for the bulk of p isotopes 
within a factor of ~3 [6, 7J. The SN shock wave induces temperatures of 2-3 GK in the 
outer (C, Ne, O) layers, sufficient for triggering the required photodisintegrations. More 
massive stellar models (M>20 M ) seem to reach the required temperatures for efficient 
photodisintegration already at the end of hydrostatic O/Ne burning [8]. The decrease 
in temperature after passage of the shock leads to a freeze-out via neutron captures and 
mainly (3 + decays, resulting in the typical p-process abundance pattern with maxima at 
92 Mo (JV=50) and 144 Sm (iV=82). 

However, the 7 process scenario suffers from a strong underproduction of the most 
abundant p isotopes, 92 ' 94 Mo and 96,98 Ru, due to lack of seed nuclei with A >90. For 
these missing abundances, alternative processes and sites have been proposed, either 
using strong neutrino fluxes in the deepest ejected layers of a SNII [yp process [9]), or 
rapid proton-captures in proton-rich, hot matter accreted on the surface of a neutron 
star [rp process jTU], Hi] ) . A few p nuclides may also be produced by neutrino-induced 
reactions during the 7-process. This "v process" [12] was additionally introduced 
because the 7 process alone strongly underproduces the odd-odd isotopes 138 La and 
180m Ta. These two isotopes could be the result of excitation by neutrino scattering on 
pre-existing s-process seed nuclei, depending on the still uncertain underlying nuclear 
physics. 

Modern, self-consistent studies of the 7-process have problems to synthesize p nuclei 
in the regions A < 124 and 150 < A < 165 [8J. It is not yet clear whether the 
observed underproductions are only due to a problem with astrophysical models or also 
with the nuclear physics input, i.e. the reaction rates used. Thus, the reduction of 
uncertainties in nuclear data is strictly necessary for a consistent understanding of the 
p process. Experimental data can improve the situation in two ways, either by directly 
replacing predictions with measured cross sections in the relevant energy range or by 
testing the reliability of predictions at other energies when the relevant energy range 
is not experimentally accessible. In this context we have carried out p-process network 
calculations with a modified reaction library which uses the most recent experimental 
and semi-empirical (ri.,7) cross sections from the "Karlsruhe Astrophysical Database 
of Nucleosynthesis in Stars" project, KADoNiS v0.2 [13]. This aims to be a step 
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towards an improved reaction library for the p process, containing more experimental 
data. However, it has to be kept in mind that the largest fraction of the p-process 
network contains proton-rich, unstable isotopes which are not accessible for cross section 
measurements with present experimental techniques. Hence there is no alternative to 
employing theoretical predictions for a large number of reactions. Typically, these come 
from Hauser-Feshbach statistical model calculations [H] performed with the codes NON- 
SMOKER [13 H6] or MOST [17]. 

2. p-process network calculations 

We studied the p process in its manifestation as a 7 process. The network calculations 
were carried out with the program "pProSim" [18]. The underlying network was 
originally based on a reaction library from Michigan State University for X-ray bursts 
which included only proton-rich isotopes up to Xenon. For p-process studies it 
was extended by merging it with a full reaction library from Basel university [19|. 
That reaction library was mainly based on NON-SMOKER predictions with only few 
experimental information for light nuclei. For the present calculations, the library was 
updated by inclusion of more than 350 experimental and semi-empirical stellar (n, 7) 
cross sections from the most recent version of the " Karlsruhe Astrophysical Database of 
Nucleosynthesis in Stars" (KADoNiS v0.2). Due to detailed balance this modification 
also affects the respective (7, n) channels. 

The abundance evolution was tracked with a parameterized reaction network, based 
on a model of a supernova type II explosion of a 25 M star [7J. Since the p-process 
layers are located far outside the collapsing core, they only experience the explosion 
shock front passing through the O/Ne burning zone and the subsequent temperature 
and density increase. Both, the seed abundances and the respective temperature and 
density profiles, were taken from [7J and are not calculated self-consistently. The p- 
process zone in the simulation was subdivided into 14 single layers. This is consistent 
with what was used in [18] and thus the results can directly be compared. 

The final 7-process abundances depend very sensitively on the choice of the initial 
seed abundance. This initial abundance is produced by in-situ modification of the stellar 
material in various nucleosynthesis processes during stellar evolution. This also means 
that the respective O/Ne layers can receive an abundance contribution from the weak 
s process (core helium and shell carbon burning during the Red Giant phase of the 
massive star) in the mass region up to ^4=90. The s-process component depends on the 
mass of the star and also on the neutron yield provided by the 22 Ne(a, n) 25 Mg neutron 
source. If the p-process layers during the explosion are located within the convective 
zones of the previous helium burning phases, the s-process abundance distribution in 
all layers can be assumed to be constant. Thus, it has to be emphasized that the 
results strongly depend on the adopted stellar model (e.g. Ref. [8j found significant 7- 
processing occurring already in late stages of hydrostatic burning). Additionally, stars 
with different masses will exhibit different s- and p-processing and true p-abundances can 
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only be derived when employing models of galactical chemical evolution. Nevertheless, 
for better comparison with previous studies, here we implemented the same approach 
as used in (6J [7J [18] (see also Fig. CP . 



3. Results and Interpretation 

The results of our simulations are shown in Fig. [1] as 'normalized overproduction factors' 
<Fi> I F [20] . This value gives the produced abundance relative to the solar abundances 
of Anders and Grevesse j2j. Ranges of variations of this factor for SN type II explosions 
with stellar masses 13 M < M star <25 M & are published in Fig. 4 of Ref. [7]. The 
factor Fn is the so-called 'averaged overproduction factor' which is a measure for the 
overall enrichment of all p isotopes. This value changed only by -0.6% from F =86.3 
(previous library [18]) to F =85.8 (modified library). The quantities <Fi> are the 
'mean overproduction factor' calculated from the mass of the respective isotope i in the 
p-process zone divided by the total mass of the p-process zone times the respective solar 
mass fraction of the isotope. In Fig. [1] our results with the modified reaction library 
are compared to the results published in [18] with the previous set of reaction rates 
to examine the influence of the newly implemented neutron capture data of KADoNiS 
v0.2. The result is a decrease in the normalized overproduction factor for almost all p 
isotopes by an average value of -7% (shown in the right part of Fig. [I] as dashed line). 
The largest deviations occur for 84 Sr (-20.5%), 136 Ce (+30.6%), 156 Dy (-39.2%), 152 Gd 
and 158 Dy (-22.3% each), 180 W (-22.9%), and 190 Pt (-32%). 
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Figure 1. Left: Normalized overproduction factors derived with the previous [18] 
(open squares) and the modified (full squares) reaction library. Additionally the values 
from a 25M0 star model of Rayet et al. [7] are given for comparison (grey stars) . A 
value equal to unity corresponds to the solar abundance. Right: Abundance difference 
for each p isotope between the modified and the previous reaction library. 



In general, these differences are relatively strong in the mass region ^4=150-170. 
Reaction flux plots of this mass region reveal that the main flux proceeds by (7,0), 
(7, n), (n, 7), and (n, a) reactions, and is lower by 1-2 orders of magnitude compared 
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to the previous reaction library. This drastic change cannot be explained by the 
implementation of larger experimental (n, 7) cross sections of some of these p isotopes, 
since this would in turn also increase the production channel via (7, n) reactions from 
heavier isotopes [20]. Surprising at first glance, the cause of the differences are changes 
in the neutron rates at higher masses. A significant fraction of the seed abundances for 
the p process is located in 209 Bi and the Pb isotopes and converted to nuclei at lower 
mass by photodisintegration sequences starting with (7, n) reactions. 

The importance of experimental data is strongly emphasized by these findings. 
Because of the magicity or near-magicity of the Pb and Bi isotopes, individual resonances 
determine the cross sections and Hauser-Feshbach theory is not applicable [2T1 122] . 
Furthermore, from discrepancies between resonance and activation measurements 
[23l El] and from theoretical considerations [25], it has been previously found that a small 
direct capture component contributes to neutron capture on 208 Pb (2T] [25] . The interplay 
of resonant and direct capture contributions is difficult to handle in theoretical models 
and experiments prove to be indispensable. This also explains why there are deviations 
up to factor of 1.7-2 between the NON-SMOKER predictions and experimental data for 
these nuclei |20j. In fact, Hauser-Feshbach models cannot be applied there. For the cases 
where the statistical model can be applied, the average NON-SMOKER uncertainty for 
stellar (n, 7) cross sections is only ±30% or even better. 

The nuclides 152 Gd and 164 Er are not produced in the present simulation based on a 
25 M Q star. This is consistent with previous work finding large s-process contributions to 
these nuclei. Also, the two odd-A isotopes 113 In and 115 Sn are not produced. The latter 
underproduction problem is known since a long time [6], [T] . The initial seed abundances 
of 113 In and 115 Sn are destroyed by the 7 process, since the destruction channel is much 
stronger than the production channel. Thus, it appears as if the nuclides 152 Gd, 164 Er, 
113 In, and 115 Sn have strong contributions from other processes and it is conceivable 
that they even may not be assigned to the group of p nuclei. 

Nemeth et al. |26j determined the contributions for 113 In, and 115 Sn with the 
(out-dated) classical s-process approach to be very small (less than 1% for both cases). 
These calculations in the Cd-In-Sn region are complicated since many isomeric states 
have to be considered, and the r-process may contribute to the abundances of 113 Cd 
and 115 In. Although these two isotopes have quasi-stable ground-states, the /3-decays 
of r-process progenitor nuclei can proceed via isomeric states: 112,9 Ag — > 113m Cd — > 
113 In and 115g Cd —> 115m In — * 115 Sn. In [26] only part of the missing abundances could 
be ascribed to post-r-process /3-decay chains, leaving rather large residues for other 
production mechanisms. In view of the progress achieved in s-process calculations using 
the TP-AGB star model, and with the availability of an improved set of (n, 7) cross 
sections it appears worth while to update these older calculations. The new reaction 
library KADoNiS [13] includes the latest Maxwellian averaged cross sections from very 
accurate time-of-flight measurements for the Cd and Sn isotopes j27J [28], [29], [30] and 
will soon be complemented by a measurement of the partial neutron capture cross 
section to the 14.1 y isomeric state in the s-process branching isotope 113 Cd with the 
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activation technique, which is underway at Forschungszentrum Karlsruhe. With this 
additional information new calculations are expected to provide more accurate s- and 
r-contributions for 113 Cd and 115 In. Based on these results, the p-contributions can be 
estimated as the residual via N p =N Q -N s -N r . 
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